Inhomogeneous isospin distribution in the reactions 
of 28 Si + 112 Sn and 124 Sn at 30 and 50 MeV/nucleon. 
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We have created quasiprojectiles of varying isospin via peripheral reactions of 28 Si + 112 Sn and 
124 Sn at 30 and 50 MeV/nucleon. The quasiprojectiles have been reconstructed from completely 
isotopically identified fragments. The difference in N/Z of the reconstructed quasiprojectiles allows 
the investigation of the disassembly as a function of the isospin of the fragmenting system. The 
isobaric yield ratio 3 H/ 3 He depends strongly on N/Z ratio of quasiprojectiles. The dependences 
of mean fragment multiplicity and mean N/Z ratio of the fragments on N/Z ratio of the quasipro- 
jectile are different for light charged particles and intermediate mass fragments. Observation of a 
different N/Z ratio of light charged particles and intermediate mass fragments is consistent with an 
inhomogeneous distribution of isospin in the fragmenting system. 
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There has been significant interest in multifragmentation of excited nuclear matter for many years. While there has 
' been some success in understanding the process of multifragmentation and describing this phenomenon in terms of a 
, liquid gas phase transition [jjj, those efforts have often treated the nucleus as a single component nuclear liquid. In fact, 
the nucleus is a two component nuclear liquid. Early work by Lamb Q laid the foundations for treating the nucleus 
as a two component system although the results were not connected to multifragmentation. Statistical calculations 
describing multifragment disassembly predicted that much of the neutron excess would be observed as free neutrons 
Thermodynamic calculations by Muller and Serot jl| lead us to the idea that, for the very neutron rich systems, 
there may exist a distribution of the excited nuclear matter into a neutron rich gas and a more symmetric liquid. 
It is predicted by both lattice gas and mean held calculations ^) that fragmentation of a system of asymmetric 
nuclear matter will express the characteristics of a two component system. This would result in a second order 
phase transition. Also within a dynamical code, constructed to study influence of charge asymmetry on spinoidal 
decomposition of nuclear matter at sub-saturation densities, differences in the fragmentation are seen as a function 
' q ■ of isospin asymmetry (?]] . 

The difference of the mean N/Z ratio of the light charged particles with Zf < 2 ( LCPs ) and of the intermediate 
mass fragments with Zf > 3 ( IMFs ) may be a possible experimental signature of a separation into a gas ( resulting 
mostly in emitted LCPs ) and a liquid ( IMFs ). An enhanced production of neutron rich H and He isotopes from 
the neck region has been seen experimentally in mid-peripheral collision s [R ^|, while a favored emission of more 
symmetric heavy clusters in the mid-rapidity region has also been shown JlQ] . Quite neutron deficient residues have 
also been seen in intermediate energy reactions |ll|, [l^]. Recent results from the reaction of 112 ' 124 Sn + 112 ' 124 Sn 
indicate that the relative abundance of free neutrons increases as the neutron content of the colliding system increases 
p3[ . The present work will demonstrate that, in a well defined system with isotopic identification of all charged 
fragments, the mean value of the N/Z ratio of LCPs is more sensitive than mean N/Z ratio of IMFs to the N/Z of 
quasiprojectile. The more asymmetric the system is the stronger it will favor breaking up into still more asymmetric 
light fragments while the N/Z ratio of heavier fragments remains relatively insensitive. 

This experiment was done with a beam of 28 Si impinging on >1 mg/cm 2 of 112 ' 124 Sn self supporting targets. 
The beam was delivered at energies 30 and 50 MeV/nucleon by the K500 superconducting cyclotron at the Cyclotron 
Institute of Texas A&M University. The detector setup is composed of an arrangement of 68 silicon - CsI(Tl) telescopes 
covering polar angles from 1.64° to 33.6° in the laboratory frame. Each element is composed of a 300 /jm surface 
barrier silicon detector followed by a 3 cm CsI(Tl) crystal. The detectors are arranged in five concentric rings. The 
geometrical efficiency is more than 90% for each ring. These detectors allow for isotopic identification of light charged 
particles and intermediate mass fragments up to a charge of Z—5. A more detailed description of the detectors and 
electronics can be found in ref. |l4j. Energy calibration of the silicon detectors was achieved with the use of a 228 Th 
source providing 6 calibration points from 5.42 to 8.78 MeV. The CsI(Tl) detectors were calibrated from the measured 
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energy loss in the silicon detector of a given particle of a certain mass and charge and from its total kinetic energy 
calculated with SRIM-96 Jig] . In addition, the punch-through points of hydrogen isotopes in CsI(Tl) crystals were 
also used. The overall procedure gives an energy resolution of typically 3%. More details of the experimental setup 
and calibration can be found in ref. ]iq] . 

The total charge of quasiprojectile ( QP ) was restricted to the values Z tot =12-15. We have required that all 
detected fragments are isotopically identified. This highly exclusive set of data possesses information on fragmentation 
of highly excited projectile-like prefragments over a range in isospin that can be measured on an event by event basis. 
Simulations have been done with a hybrid code treating the two stages of the reaction - nucleon exchange and 
quasiprojectile deexcitation. The first stage of the calculation was simulated with the Monte Carlo code of Tassan- 
Got p7j |. This code implements the model of deep inelastic transfer jOt. For the deexcitation of the highly excited 
quasiprojectile we used the statistical model of multifragmentation ( SMM ) |l9| . The quasiprojectile event sequences 
generated by the DIT code have been used as input of the SMM simulations. The results of the calculation were 
filtered through the FAUST software replica. This hybrid calculation reproduces very well the quasiprojectile and 
fragmentation observables in these reactions [^o) . 

In figure [j] we show the effect of including emitted neutrons in the calculation of N/Zqp. Since the experimental 
setup cannot detect free neutrons it is important to see if this might bias the data. The N/Z„ onetlt is the neutron to 
proton ratio of the quasiprojectiles simulated using hybrid DIT/SMM simulation and filtered using FAUST software 
replica. Thus N/Z„ oneut applies to quasiprojectiles where only charged particles are included. The N/Z neut is the back- 
traced total neutron to proton ratio of the primary quasiprojectiles for the same set of events. The line represents equal 
values for N/Z„ etlt and N/Z noneut . The data presented on figure |l| represent the reaction 28 Si(50MeV/nucleon)+ 112 Sn. 
While there is some deviation from the line due to neutron evaporation, the calculation predicts that N/Z„ eut tracks 
very well with N/Z norieut . The situation is similar also when using 124 Sn target or lower projectile energy. Therefore 
for the remainder of this paper N/Zqp will be N/Z n0 neut - 

In figure|| ^ s snown the dependence of the isobaric yield ratio 3 H/ 3 He on N/Zqp. The yield ratio 3 H/ 3 He increases 
as N/Zqp increases. Similar effects have been seen in other systems where isotopic or isobaric ratios have been 
plotted as a function of N/Z of the reacting system || ^l). In the present data, the N/Zqp is the reconstructed N/Z 
of all detected charged particles. If one fits the data with an exponential function exp(ai + 0,2 N/Zq P ), the slope 
0,2 is nearly identical for the reaction with the 112 Sn target as with the 124 Sn target. The slopes presented in table 
I show that this increase in the yield ratio 3 H/ 3 He as a function of N/Zqp is much larger than would be expected 
from the change in N/Zqp alone. While the integrated values on the two targets are different when broken down by 
the N/Zqp the behavior is similar. If one looks at the difference with bombarding energy, the dependence of yield 
ratio 3 H/ 3 He on N/Zqp is lessened at the higher energy. This is consistent with the lattice gas calculation that 
predicted that the 3 H/ 3 He ratio would be enhanced at lower temperatures for the system with excess neutrons. At 
relatively high temperatures above 12 MeV, the ratio would asymptotically approach the N/Zqp value. The hybrid 
DIT/SMM calculations reproduce the experimental yield ratio 3 H/ 3 He rather well. This shows that statistical model 
of multifragmentation correctly describes the influence of isospin on the production rates of fragments. The overall 
3 H/ 3 He ratio for the reaction with the 112 Sn target is 0.25 and with the 124 Sn target is 0.35 at 30 MeV/nucleon and 
0.44 and 0.62 respectively at 50 MeV/nucleon. This increase in the yield of neutron rich fragments with an increase in 
the N/Z of the reacting system is consistent with other studies of isobaric ratios [|l3, 21 . What figure || demonstrates 



is that this is an integrated feature and merely shows that the system created with the neutron rich target is more 
neutron rich. But when one takes the 3 H/ 3 He ratio as a function of the N/Zqp, the behavior may be seen much more 
clearly. 

In figure || we have investigated the effect of excitation energy on the 3 H/ 3 He ratio. The solid squares represent 
all events, where the excitation energy is less than the mean excitation energy per nucleon of the events included 
in the figure ||c ( 28 Si+ 112 Sn at 50 MeV/nucleon, e* < 5.5 MeV/nucleon ). The open squares represent the events 
whose excitation energy is greater than the mean excitation energy. The slope is steeper for events at lower excitation 
energies. We observed a similar behavior for all reactions regardless of the bombarding system or beam energy. The 
results of hybrid calculation ( lines ) are again in reasonable agreement with experiment. Further work investigating 
the excitation energy dependence of 3 H/ 3 He ratio is forthcoming p2| . 

In figure [|a we show the multiplicity of charged particles versus N/Zqp, again for the reaction 28 Si+ 112 Sn at 50 
MeV/nucleon. The squares represent the multiplicity of all charged particles. The multiplicity of charged particles 
increases as N/Zqp decreases. This is then broken down into the multiplicity of light charged particles ( circles ) and 
the multiplicity of intermediate mass fragments ( triangles ) . Here we can see that the multiplicity of LCPs increases 
with decreasing N/Zqp. Meanwhile the multiplicity of IMFs increases as N/Zqp increases. The situation is similar 
for all cases. The results are identical for different targets at the same projectile energy. The dependence of the 
multiplicity of IMFs on N/Zqp is practically the same for all cases. The increase in overall multiplicity of charged 
particles with increasing beam energy is similar to the increase of the multiplicity of LCPs. The results of the hybrid 
calculation ( lines ) are consistent with the experiment. Again the presented information is consistent with, but much 
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more illustrative than, previous works where the difference in isospin of the excited system was only approximately 
known on average Q . Previous works show that the multiplicity of IMFs as a function of multiplicity of charged 
particles increases for the more neutron rich system. While this is true, the present data shows that this is concurrent 
with the decrease in the multiplicity of LCPs for neutron rich systems. Our data shows that this effect weakens at 
higher energies, possibly towards the disappearance, that was seen by Miller et al. |24j] at much a higher energy. Our 
data would be consistent also with the temperature dependence predicted by lattice-gas model calculations [|| . 

In figure |]b we show the ratio jy/f^p as a function of N/Zqp. The ratio jv/z^p decreases with increasing N/Zqp. 
As there are fewer neutrons available the excess protons go into the smaller fragments rather than the larger fragments. 
The least neutron rich quasiprojectiles with N/Zqp w 0.5 prefer to breakup into very neutron deficient LCPs and much 
more symmetric IMFs. This inhomogeneous isospin distribution washes out rapidly as N/Zqp approaches the region 
of /3-stability. Such a behavior may be understood when taking into account that the N/Z range of detected LCPs 
is much wider than the N/Z range of IMFs where usually only few stable and nearly stable isotopes are produced. 
The observed inhomogeneous isospin distribution is likely caused by more favorable energy balance of the deposition 
of proton excess into LCPs, either free protons or light proton rich clusters ( see Figs. §,[f )• When extrapolating this 
trend towards very neutron rich quasiprojectiles the isospin distribution would most likely reappear again in the form 
of neutron rich LCPs and more symmetric IMFs. Such an extrapolation is consistent with the predicted asymmetric 
liquid gas phase transition [Q where in addition to the inhomogeneous isospin distribution a spatial separation of 
fractions with different isospin occurs. For neutron rich nuclei, such a phase transition is predicted at the values of 
N/Z exceeding 1.5, corresponding to proton concentrations 40 % or less. The data presented here raises the question 
whether an analogous separation of the system into subsystems with different isospin may be expected for very neutron 
deficient fragmenting systems. Our data support such a scenario. Nevertheless, one has to take into account that the 
studied fragmenting system is quite small and the number of possible exit channels may play an important role in 
the determination of final partition of fragments. The influence of the secondary decay was estimated by simulations 
where only an emission of hot fragments took place. The number of simulated events, where all fragments have 
been isotopically resolved, decreased for a given overall number of simulated events by 50 % at projectile energy 30 
MeV/nucleon and by 90 % at 50 MeV/nucleon. This may be explained by a lower mean multiplicity and higher mean 
mass and charge of hot fragments, what in our case leads to a lower number of simulated events where only fragments 
up to boron ( Zf < 5 ) are emitted. Nevertheless, the absence of secondary emission stage did not influence the ratio 
n/z IM f ag a f unc tion of the N/Zqp, which remained unchanged within statistical errors. This implies that the isospin 
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dependence of the ratio jy/zlcp is, according to the statistical model of multifragmentation, determined in the hot 
multifragmentation stage. The neutron emission does not seem to influence the presented data dramatically, as may 
be seen in figure [l| 

We have demonstrated that systems of varying isospin can be created and studied by projectile fragmentation 
reactions. The spread in N/Z of the fragmenting system is larger than the difference between the reacting systems. 
This can be observed by the complete isotopic reconstruction of the quasiprojectile. We have investigated the isobaric 
yield ratio 3 H/ 3 He as a function of N/Z ratio of the fragmenting system. The yield ratio 3 H/ 3 He increases as N/Zqp 
increases. This depends more significantly on the fragmenting system than on the reacting system. At higher energies 
this dependence is lessened. We conclusively demonstrated that the multiplicity of charged particles depends on the 
isospin of the fragmenting system. As the N/Zqp decreases the overall multiplicity of charged particles increases. 
The multiplicity of LCPs increases dramatically while the multiplicity of IMFs decreases. This effect is less significant 
at higher bombarding energy. We have also presented evidence of an inhomogeneous isospin distribution of the 
fragmenting system into two fractions with different values of N/Z. 
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FIG. 1: Average N/Z ratio of all fragments emitted from the quasiprojectile versus those bound in charged fragments for 
the reaction 28 Si(50MeV/nucleon)+ 112 Sn. Symbols represent the hybrid DIT/SMM calculation. Solid line indicates the case 



FIG. 2: Isobaric yield ratio 3 H/ 3 He versus N/Zqp for the reactions of 28 Si beam with tin targets, (a) - 30MeV/nucleon, 112 Sn, 
(b) - 30MeV/nucleon, 124 Sn, (c) - 50MeV/nucleon, 112 Sn, (d) - 50MeV/nucleon, 124 Sn. Symbols represent the experimental 
data. Lines represent the hybrid DIT/SMM calculation. 

FIG. 3: Isobaric yield ratio 3 H/ 3 He versus N/Zqp for the reaction 28 Si(50MeV/nucleon)+ 112 Sn. Symbols represent the 
experimental data for two different cuts on apparent excitation energy per nucleon ( Solid squares - e* < 4eam open squares 
- e* > 4e«n )• Lines represent the hybrid DIT/SMM calculation. ( Solid line - e* < 4eam dashed line - e* > e* mean ). 

FIG. 4: (a) - Multiplicity of charged fragments (squares), LCPs (circles) and IMFs (triangles) versus N/Zqp. Corresponding 
lines represent the hybrid calculation, (b) - Experimental ratio of the mean values of N/Z of LCPs and IMFs versus the N/Zqp. 
Data are given for the reaction 28 Si(50MeV/nucleon)+ 112 Sn. 



TABLE I: Values of the slope parameter obtained from the fit log 10 ( , 3 tt , ) = ai + o,2{^)qp- 



Reaction 


Projectile energy 


ci2(exp) 


a2(calc) 




[MeV/nucleon] 






28 Si+ 112 Sn 


30 


3.99±0.39 


3.89±0.16 


28 Si+ 124 Sn 


30 


4.13±0.11 


4.57±0.24 


28 Si+ 112 Sn 


50 


3.16±0.18 


3.15±0.16 


28 Si+ 124 Sn 


50 


2.83±0.27 


3.10±0.11 



